Generalized Effective Potential Landau Theory for Bosonic Quadratic Superlattices 
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We study the properties of the Bose-Hubbard model for a quadratic optical superlattice. To this 
end we generahze a recently established effective potential Landau theory for a single component to 
the case of multi components and find not only the characteristic incompressible solid phases with 
fractional filling, but also obtain the underlying quantum phase diagram in the whole parameter 
region at zero temperature. Comparing our analytic results with corresponding ones from quantum 
Monte Carlo simulations demonstrates the high accuracy of the generalized effective potential Lan- 
dau theory (GEPLT). Finally, we comment on the advantages and disadvantages of the GEPLT in 
view of a direct comparison with a corresponding decoupled mean-field theory. 
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I. INTRODUCTION 



Systems of ultracold bosonic gases in optical lattices 
have recently become a major field in physics researchfi"— 
After their theoretical suggestion^ and first experimen- 
tal realization using counter-propagating laser beams^ it 
soon became clear that they establish a versatile bridge 
between the field of ultracold quantum matter and cor- 
related condensed matter systemsi^ 

One of the most famous examples is the Bose-Hubbard 
model^i^ which undergoes a quantum phase transition 
from a Mott insulator to a superfluid phase due to the 
competition between the atom-atom on-site interaction 
and the hopping amplitude. This transition can be 
demonstrated experimentally by time-of-flight absorp- 
tion pictures?^ or measuring the collective excitation 
spectra via Bragg spectroscopyi^i^ Recent research ef- 
forts have targeted more complex systems, which include 
long range interactions (e.g. from dipolar bosons^Siii) , 
mixtures of several componentsi^ii^ and more interest- 
ing lattice geometries, such as frustrating or superlattice 
structures ji^"— Accordingly, the corresponding phase di- 
agrams become richer and more complex, which includes 
the possibility of phases with crystalline order or super- 
solidity. A crystalline density wave phase generally oc- 
curs at fractional filling and it has been proposed that 
the corresponding commensurate order could be detected 
by measuring correlations with time-of-flight and noise- 
correlation techniques.— Furthermore, recent experimen- 
tal progress in achieving single-site addressability in op- 
tical lattice structuresi^Ti^^ nourishes the prospect to di- 
rectly observe the density waves in the near future. 

One of the simplest ways to induce density wave order 
is to create a superlattice with additional lasers r^"— so 
that the potential depth is slightly different on one sub- 
lattice, while the interaction strength and hopping am- 
plitude will remain almost uniform. The corresponding 
Bose-Hubbard model Hamiltonian on a square or cubic 



lattice is given by2& 
U 
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where A/i stands for a small additional chemical potential 
on sublattice A compared to sublattice B, as illustrated in 
Fig.[T] As we will show this model exhibits an interesting 
competition between Mott and density wave phases. 

From a theoretical point of view the study of in- 
teracting bosons and quantum phase transitions is far 
from trivial)^ The possible phases in different kinds 
of optical superlattices have so far been analyzed by 
numerical approaches^"— and decoupled mean field 
theory,— i^i^i^^ which are able to predict the phase di- 
agrams qualitatively. However, it is known that mean 
field theory or variational calculations can have signifi- 
cant deviations from the unbiased high-precision numer- 
ical results 1^ The purpose of this paper is to present 
a reliable quantitative method to determine non-trivial 
phases of multi-component boson systems based on re- 
cent advances to use a systematic Landau theory with 
an effective potential that can be estimated quantita- 
tively from the microscopic model, e.g. by diagram- 
matic methods.—"— Whereas the first hopping order 
of the effective potential Landau theory leads to sim- 
ilar results as mean-field theory^i higher hopping or- 
ders have recently been evaluated via the process-chain 
approach)^"— which determines the location of the 
quantum phase transition for the single component Hub- 
bard model for cubic as well as triangular and hexag- 
onal optical lattices to a similar precision as demand- 
ing quantum Monte Carlo simulations4^i^ Thus, it be- 
comes even possible to calculate the critical exponents 
of the corresponding quantum phase transition4ii^ We 
now present a generalized effective potential Landau the- 
ory (GEPLT), which extends those concepts to multi- 
component systems and to phases with non-trivial crys- 




FIG. 1: (Color online) Schematic illustration for a quadratic 
superlattice in two dimensions. The solid line square repre- 
sents one type of the unit cell. The solid curve (green) shows 
the optical potential in the x direction. Lattice sites A are 
deeper by A/i than lattice sites B. 



talline order parameters. In particular, for the model in 
Eq. ([T]) the GEPLT approach gives excellent quantitative 
estimates for the location of the phase boundaries com- 
pared to unbiased quantum Monte Carlo simulations. 

At first, we briefly review the effective potential Lan- 
dau theory for the single component Bose-Hubbard 
model in Sec. |lTl Then, we extend this method step by 
step from one component to the general superlattice case 
in Section Uni After that, we apply this GEPLT method 
to the simple superlattice model in Eq. ([T|) and determine 
the resulting quantum phase diagram at zero tempera- 
ture in the whole parameter region in Section IIVI Both 
the advantages and disadvantages of GEPLT are revealed 
by comparing it with a decoupled mean-field theory in 
Section |Vl Finally, Section I VII provides the conclusions 
and sketches related problems in an outlook. 



II. EFFECTIVE POTENTIAL LANDAU 
THEORY 

Let us first consider the Bose Hubbard model in Eq. ([T]) 
for the well-studied case of A^ = OA^ The second-order 
quantum phase transition between the Mott insulator, 
which occurs for t/U <^ 1, and the superfluid, which 
is realized for t/U 3> 1, is intimately connected with a 
spontaneous breaking of the underlying U(l)-symmetry 
of the Bose-Hubbard model ([l} . To describe this theoret- 
ically, we transfer the usual field-theoretic approach for 
thermal phase transitions^Si^ to quantum phase transi- 
tions and couple the creation and annihilation operators 
to external source fields with uniform strength J and J* 
within a Landau theor y'^^i^^ 



^bh( J, J*) = Hbh + J2 (-^"l + •^*«. 



(2) 



The transition from the Mott insulator to the superfluid 
phase is described by the emergence of a non- vanishing 



order parameter which is defined due to homogeneity ac- 
cording to ?/> = (fli), ip* = (al). The free energy corre- 
sponding to ^ 



F{J, J* 
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allows to determine this order parameter via 
1 dF{J,J*) ,,_ 1 dF{J,J*) 



V' 



Ns dJ* 



r = 



N, 



dJ 



(3) 



(4) 



where Ns denotes the number of lattice sites. Equation 
(|4]) motivates that it is possible to formally perform a 
Legendre transformation from the free energy F{J,J*) 
in order to arrive at an effective potential T^^jip*) that 
is useful in a quantitative Landau theory 



T{'^,^P*)=F/Ns~ib*J-ijJ\ 



(5) 



Due to Legendre identities the external sources can be 
reobtained from derivatives of the effective potential 






-J, 






-J* 



(6) 



The original Bose-Hubbard Hamiltonian ([T]) is restored 
from ([2]) for vanishing currents, i.e. by setting J = J* = 
0. In this limit we conclude from ([5|) that the effective po- 
tential reduces to the free energy. Furthermore, Eq. 1^ 
then implies that the order parameter tp, tjj* of the system 
follows from extremizing the effective potential. A triv- 
ial extremum t/i = corresponds to the Mott-insulator 
phase, whereas a non-vanishing extremum "0 7^ occurs 
in the superfluid phase. 

The free energy ^ reduces at zero temperature to the 
ground-state energy, which can be calculated in a power 
series of both the hopping parameter t and the source 
terms J, J* by using the Rayleigh-Schrodinger perturba- 
tion theoryj^^— Due to the underlying U(l)-symmctry 
of the Bose-Hubbard Hamiltonian ([ij the expansion is 
only a power series in terms of | Jp 



F{J,J\t) = Ns Foit) + J2c2pit) \J 



2p 



(7) 



where the respective expansion coefficients are accessible 
via a hopping expansion 



C2p(0 = E(-^)'"": 



2n„(n) 
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From Eqs. (H]), ([5]), and ([7]) we then obtain the effective 
potential of the Bose-Hubbard Hamiltonian ([T]) in the 
following pcrturbative form 



Ti^,r,t) = Fo{t) 



1 



■\M' + ^\i'\' 
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According to the Landau theory for second-order phase 
transitions, the critical line between the Mott insulator 



and the supcrfluid phase follows from finding the zero of 
the second-order coefficient in ^. In order to solve the 
resulting equation l/c2(tc) = 0, we expand it in a power 
series of the hopping parameter t 
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Thus this gives us an algebraic equation for tc^ whose 
degree depends on the respective hopping order which 
is taken into account. The number of the roots is the 
same as the order of i, but only the smallest real positive 
root is identified as an appropriate approximation for the 
location of the quantum phase transition. 

As mentioned in the introduction, the effective poten- 
tial Landau theory was quite successful in calculating 
the quantum phase boundary for the single component 
system<^"— However, it cannot be used to treat a su- 
perlattice system, since more than one order parameter 
appears. Therefore, we will work out in the next section a 
corresponding extension to multi components which over- 
comes this problem. 



III. GENERALIZED EFFECTIVE POTENTIAL 
LANDAU THEORY 



In order to describe a superlattice or a multi- 
component system, we have to introduce several sites or 
degrees of freedom at each lattice point. In other words, 
we introduce a larger unit cell at each lattice point, la- 
beled by j, together with a basis of size m, labeled by 
I = l,..,m. The generalized Bose- Hubbard Hamiltonian 
with m bosons in each unit cell is therefore given by 

^BH = "X! H [*j(/)j'(i')«l(o«j'(;')+h-c.] (11) 

where ajin denotes the boson annihilation operator at 
lattice point j with basis index (/). Hopping ij(;).j'(;') 
can occur between any basis and lattice position, while 
the repulsion Vn) s-cts for now only between bosons of 
the same lattice point and basis index. The chemical po- 
tential /^(i) depends on the basis index which is analogous 
to sublattices A and B in Eq. (jlj. 

We now model the symmetry-breaking by introducing 
the source vectors J = (Ji, . . . , JniS^ ,1^ = {Ji, ■ • ■ , Jm) 
according to 



H, 



BHl 



H, 



BH 
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• (12) 



By generalizing the procedure from a single component to 
multi components, we use perturbation theory in order to 
determine the free energy at zero temperature in a power 
series of both the hopping parameters ij(i),j'{V) and the 
source vectors J, J^. In principle, we need an expansion in 
terms of all relevant hopping parameters tj{i),j'(i')^ but to 
illustrate the process we consider here the case that only 
one hopping element t dominates (e.g. between nearest 
neighbors) and all others are neglected 



F{Sj\t)^N, Fo{t)+PC2{t)J + - 



(13) 



The matrix elements c^n, (t) of C2 (i) are then given by a 
hopping expansion of the form 



=2«'W-EH)""2«^ 



(14) 



n=0 



The order parameter vectors give different values for each 
basis index, but are independent of lattice points ^ = 
(V'i,...,V'm)'^,^'^ = (V'r,---,V'm) according to 



* 

^t = 



((«i 



(flm)) 



^l/i • • • I \^mi 



and we observe 



1 dF . _ 1 dF 



(15) 



(16) 



Again this motivates to perform the Legendre trans- 
formation of the free energy. The generalized effective 
potential then depends on the order parameter vectors 



r(*,*t^t) = F/Ns -jj^-* 



(17) 



Legendre identities allow to write the external sources as 
derivatives of the effective potential 



dr 

9* 



dr 



(18) 



so the order parameter vector is determined by extremiz- 
ing F in the physical limit that the external source vectors 
J, JJi^ vanish. 

Due to Eqs. (HH), dUl), and ([TTI) the effective potential 
of the system is of the form 



r(*, ^\t)^ Fait) + ¥A2it)^ + ■■■ 



(19) 



The resulting relation between the matrices A2 and C2 
can be deduced in following way. By inserting p9p into 
dUl), we get 



(20) 



Combining this with (jTS)) and (|T6)) . we read off 



(21) 



Thus, the matrix A2 turns out to be the inverse of — C2. 
As all matrix elements of C2 are given by a hopping ex- 
pansion of the form (fT4|) . we get a corresponding hopping 
expansion for each element of A2 . In matrix form the first 
terms of this hopping expansion read 
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(22) 



which reduces for a single component to pU)) . The crit- 
ical line, where the order parameter vector ^ changes 
from zero to non-zero, follows then from extremizing the 
effective potential ([T9l). In case that all components of 
the order parameter vector ^ arc non-zero, we obtain 



Det A2 = . 



(23) 



But it could also happen that only a subset of compo- 
nents of the order parameter vector VP is non- vanishing, 
which yields the condition that the determinant of the 
corresponding submatrix of A2 vanishes. The physically 
realized quantum phase boundary corresponds then to 
the smallest value of the hopping parameter t, which fol- 
lows from all these conditions. In the next section we will 
study along these lines the most simple case of a super- 
lattice system which is provided by bosons on a quadratic 
or cubic superlattice given in Eq. ^. 



IV. QUADRATIC OPTICAL SUPERLATTICE 

Similar to the continuous translational symmetry 
breaking artificially introduced by the optical lattice to 
mimic a real crystal, the optical superlattice can break 
the discrete translational symmetry to study the multi- 
components system. Beside that, it also can be used as a 
platform for disorder— and topological orde r^^i^^ prob- 
lems. Here, we apply the GEPLT to the simple quadratic 
case. 



where the symmetry C2AB = C2ba holds. Then, after 
the Legendre transformation P7)) . we obtain the effective 
potential p^ . where we have ^ = {{o-a), {'Ib})'^ and A2 
is the inverse of —C2 according to Eq. (PT|) . When the 
second-order quantum phase transition occurs, the van- 
ishing order parameter vector ^ = (0, 0)""" changes from 
stable to unstable. It turns out that the smallest critical 
hopping parameter results from the condition (|23p that 
the determinant of A2 vanishes. With this we obtain in 
second order of t the following equation for the location 
of the quantum phase boundary 



(/3(0))2_ ^(0)^(1)^ 



where the abbreviations /S*^"' 



/3(2)i2 
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have been introduced. Taking into account the smallest 
root then yields 



/3(o) r ^(1) + ^(/3(i))2 + 2/3(2) 



/3(2) 



(27) 



Thus, the problem of finding the quantum phase bound- 
ary has been reduced to the calculation of the pertur- 
bative coefficients a^l, in the respective hopping order. 
According to Appendix |A] this perturbative calculation 
can be systematically performed by using a suitable dia- 
grammatic representation. We use the unperturbed en- 



ergies 



£;(0)(nA,nB) 



— UA (J^A - 1) + -^"-s [nB ~ 1) 

- (/x -f A^) UA - [iUB , (28) 



A. Application of Effective Potential Theory 

Following the GEPLT from the previous section for the 
model in Eq. ^ we need to use for the two sublattices 
two independent source terms JJ = (Ja, t/s)"^, yielding 



H^ 



BHl 






Y.{Jb 



Jga 
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The free energy of the system can then be written as p3l) , 
where C2(i) represents a 2x2-matrix with the following 
hopping expansion 



C2AA C2AB 
C2BA C2BB 
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to define the energy differences between different particle 
number sectors 
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±(p) 
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X±W ^ [E^o) ^nA,nB) - i?(°) {uA, ub ± p)] /N,. (29) 



For p 



±1 the short notation A^ 
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^A ' ^B 
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is used. In zeroth and first hopping order we obtain the 
following results for the respective coefficients a. 
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FIG. 2: (Color online) The quantum phase diagram of a 
bosonic 2d quadratic superlattice in the whole parameter re- 
gion from first hopping order GEPLT, and the phase diagram 
projected in the t-fj, plane comparing the quantum Monte 
Carlo simulation (dotted line), the 1st order (dashed line) 
and the 2nd order (solid line) analytic results at Afi/U = 0.5. 



whereas in second hopping order we get 






z{z- l){a2AA?(ABB + Z 
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and analogous for cC^bb '^ith the indices A and B inter- 
changed. 



B. Quantum Phase Diagram 

In order to got the whole quantum phase diagram, we 
study at first the i = contribution of the effective po- 
tential in Eq. (HH), i.e. i^o(i = 0) = NsE'^^'i {nA,nB) with 
Eq. (pS)) . We assume A/i to be in the region of [0,C/). 
Similar to the normal Bose-Hubbard model, there exist 
Mott insulator phases (Mott-n), which are characterized 
by the uniform filling ha = "-b = n. However, due to the 
local offset A/i, this happens only in the regions 
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FIG. 3: (Color online) The quantum phase boundaries of a 
bosonic 3d cubic superlattice at A^/U — 0.5 which is ob- 
tained by 1st order (dashed line) and 2nd order (solid line) 
generalized effective Landau potential theory, and the quan- 
tum Monte Carlo simulation (dotted line) in the thermody- 
namic limit. 



On the other hand the density wave phases (DW-n) break 
the translational order as they have the property ns = 
ua — 1, nA = n, yielding the filling factor n + 1/2, and 
minimize the free energy in the other regions 



DW-n: {n-l)U -AiJ< ^i<{n-l)U. 



(34) 



Mott-n: U {n - I) < n < Un - A^i . 



(33) 



Hence, depending on the chemical potential offset A/i, 
we find a natural competition between Mott phases and 
density wave phases. 

Turning on the hopping processes, the quantum fluc- 
tuations will melt the different insulating phases, and the 
critical lines are determined in second hopping order by 
Eq. (P7)) after substituting the respective strong-coupling 

coefficients 02"// fromEqs. (pni) - ([5^ . The resulting quan- 
tum phase diagram for quadratic and cubic superlattices 
are shown in Fig. [2] and Fig. [3l respectively. 

From the GEPLT calculation we find for the special 
case Afi = that the Mott-1 lobe coincides with the 
single component method from Ref. [35| as expected. In 
addition, when A/i is larger than zero, the DW phase 
appears, and its region increases with increasing A/i, 
whereas the region of the Mott phase decreases corre- 
spondingly. This is a direct consequence of the transla- 
tional symmetry breaking due to the superlattice struc- 
ture. Furthermore, this observation is confirmed by a di- 
rect comparison of the lobe maxima according to Fig. |4l 
where the tips of the Mott lobes decrease with increasing 
A/t, whereas the tips of the DW lobes increase. 

Comparing the quantum phase diagram in different 
dimensions, we notice that not only the lobes of the 
Mott insulators but also the DW phases are smaller 
in three than in two dimensions, which indicates that 
the dimensionality has a similar effect on both incom- 
pressible phases. In addition, we find that the results 
of the GEPLT match well with high-precision quan- 
tum Monte Garlo simulations^"—. This indicates the 
efficiency and high accuracy of our algorithm. Note 



0.04 

0.035 

0.03 

0.025 

5 0.02 

0.015 

0.01 

0.005 

0, 





— Mott-l 

Mott-ll 

-■-Mott-lll 

— DW-I 
---DW-II 
--DW-III 

DW-IV 


-—p-^ 








/^^ ;';:::::--'-^-':::'l 




^ ,-''' ^~_~-~-"^ "■■-"■-' — 


A-'- 


■:'-'■"' 









0.1 



0.2 0.3 

An/U 



0.4 



0.5 



FIG. 4: (Color online) Maxima of Mott-n and DW-n lobes as 
a function of A/i for a 3d cubic superlattice. 



that the Bose-Hubbard model in the superlattice system 
can also be analyzed by using the decoupled mean-field 
theoryi^i^i^ii^j whose quantum phase boundary turns 
out to agree with our first-order hopping result. There- 
fore, we compare in the next section the advantages and 
disadvantages of GEPLT with this decoupled mean-field 
theory. 



V. COMPARISON WITH DECOUPLED 
MEAN-FIELD THEORY 

In order to treat a superlattice system with the de- 
coupled mean- field theory, the operators c^iiai) are de- 
composed into the mean fields ?/;* ("01), which are identi- 
fied with the order parameters, and the remaining op- 
erators i5aj((5ai), which describe the quantum fluctua- 
tions around the mean fields. Then, after neglecting sec- 
ond order terms of the quantum fluctuations and assum- 
ing due to homogeneity that the order parameters arc 
equal in the same subsystem, the Bose-Hubbard Hamil- 
tonian ([1]) decouples into a mean-field Hamiltonian on 
two subsystemsiSi^^i^i^: 



-a\^B - a^B^A - aA'ip*B - asV'AJ i (35) 



with Ho from (jA3p . Because the order parameters are 
tiny near the boundary of the second-order phase tran- 
sition, the corresponding free energy can be Taylor ex- 
panded with respect to both order parameters 
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MF 






+ f2B\lpB\ + f2ABlpAi>*B 



(36) 



where the leading term /o is equal to the leading term 
Fo{t) of GEPLT at t = 0. Thus, from the previous anal- 
ysis on Fo{t = 0) = iV,,£:(°) {nA,nB) with Eq. (HHD, we 
obtain the restrictions (P5)) and ([M|) for the chemical po- 
tential in the Mott insulator and density wave phases, 



respectively. As we only consider the system at zero tem- 
perature, the free energy is equivalent to the ground-state 
energy, which can be calculated perturbatively in the oc- 
cupation number representation. With this we get the 
second-order coefficients 
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With the conditions (j33]) and (|M)) we find for both 
second-order derivatives the inequalities 



d^F 
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f2A<Q, 



= hB < . 



(38) 



This contradicts with the minimum condition which re- 
quires that both second-order derivatives are positive at 
ipA = ipB = 0. We consider this to be a general problem 
of the multi-component decoupled mean-field theory, be- 
cause it also happens in other systems such as Kagome 
and triangular systems. Note that it can be shown that 
a single-component mean-field theory does not have this 
minimum problem. 

In order to proof that the GEPLT docs not suffer from 
such a problem, we conclude at first from Eq. (j25p 



a2AA a2AB 
a2BA a,2BB 
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Considering, for instance, the first-order result, we have 

''^"2AA"2BB ^ ^' ^2AA — 



in the Mott lobe Cj^^ 



"2AA < 0' 
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^2BB 



'^2BB ^ ^^ SO ^^ §®* ^P ^'^ ^^'S* '^^'~ 

der ^C2AAC2BB - C2AB > 0. Considering ^/c^aactbb + 
C2AB = is the phase boundary and t is decreasing from 
the supcrfluid phase to the insulator phase, the denom- 
inator of Eq. pop is positive in the insulator lobe which 
means 



and 
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(41) 



For the same reason, a2BB is also positive. Thus, the 
effective potential is really a local minimal at the zero 
point. Thus, in comparison with the decoupled mean- 
field approach GEPLT has the decisive advantage to be 
consistent for superlattice systems. 



Another advantage of our method is its higher aecu- 
racy. In comparison with quantum Monte Carlo simula- 
tions, the error of the GEPLT is less than 3% in second 
hopping order. And, according to our knowledge, such 
high accuracy is hard to reach by using other analytic 
methods. It can only be surpassed by higher hopping 
orders which could be evaluated via the process-chain 
approach of Refs. f40l - |4i |. 

However, the GEPLT also has its disadvantages. For 
the two-dimensional quadratic superlattice system we 
can not get the full lobe of the phase boundary of the 
DW-I phase in the parameter range A/x < 0.35[/ in sec- 
ond hopping order, because the radicant of the square 
root in the phase boundary Eq. (P?]) becomes negative in 
second hopping order. We suspect that this is an artifact 
of truncating the hopping expansion at second order and 
expect that this could be corrected by obtaining higher 
hopping orders. 



VI. CONCLUSION AND OUTLOOK 

In this paper, we extended the single-component ef- 
fective potential Landau theory to the general case of 
a multi-component GEPLT method. In order to in- 
clude several order parameters, we introduced the source 
vectors into the general multi-component Bose-Hubbard 
Hamiltonian. After performing the Legendre transfor- 
mation of the free energy, we obtained a generalized ef- 
fective potential, which can be determined in an expan- 
sion in hopping matrix elements. This method can be 
applied to the bosonic quadratic and cubic superlattice 
systems yielding high accuracy results for the phase dia- 
grams in second order hopping compared to QMC sim- 
ulations. Apart from the Mott insulator phases, we also 
found competing DW phases with fractional filling fac- 
tors which are induced by the translational symmetry 
breaking of the superlattice system. The dimensionality 
has a similar effect on the Mott insulator and the DW 
phases. Compared with the decoupled mean-field theory, 
the GEPLT has a higher accuracy and does not suffer 
from the local minimum problem. However, GEPLT also 
has a problem in calculating the whole quantum phase 
diagram for the DW-I phase which should be solved by 
considering higher order hopping corrections. 

As the GEPLT turned out to be a general method for 
detecting second-order quantum phase transitions in a 
system with multi order parameters, we think it is also 
suitable for frustrated superlattice systems, such as the 
triangular and the Kagome lattice. Since the supersolid- 
solid phase transition for hard-core bosons is found to be 
of second order in the triangular lattice^, the GEPLT 
introduces a promising way to detect the quantum phase 
transition in both positive and negative hopping process 
regions. Furthermore, extending this work for finite tem- 
peratures and investigating the universal properties near 
the quantum phase boundary, are certainly worth for 
more detailed studies in the future. 
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Appendix A: Strong-Coupling Peturbation Theory 



The pcrturbative coeflficients a^^, follow at zero 
temperature from applying Rayleigh-Schrodinger 
perturbation theory using a suitable diagrammatic 
representation.— By denoting the creation (annihila- 
tion) operator with an arrow line pointing into (out of) 

the site, each pcrturbative contribution of 02"// can be 
sketched as an arrow-line diagram which is composed 
of n oriented internal lines connecting the vertices and 
two external arrow lines. The vertices in the diagram 
correspond to the respective lattice sites, oriented 
internal lines stand for the hopping process between 
sites, and the two external arrow lines are representing 
creation and annihilation operators, respectively. Table |T] 
presents all non-vanishing arrow-line diagrams as well 

(n) 

as the associated multiplicities aj//' ^P to the second 
hopping order. 

Note that the arrow-line diagrams only depict the pos- 
sible hopping processes. In order to determine each non- 
zero pcrturbative term a^^^ii, we also invoke a line-dot 
diagrammatic representation which has been worked out 
for a single component method in Ref. [S^- To this end 
we consider the general situation that a Hamiltonian H 
decomposes into an unperturbed term Hq and two per- 
turbation terms V, W , i.e. 



H = Ho + XV + SW. 



(Al) 



where A and 5 are small parameters. We then calculate 
the zero-temperature free energy by using perturbation 
theory. Each term is related to several line-dot diagrams 
which stem from the following rules: 

• The dots labeled by 1 and 2 represent the pcrtur- 
bative terms V and W, respectively. 

• The internal lines connecting two adjacent dots are 
associated with the factor 






1 



-E, 



(0) 



_|^(0)\/^(0)| 



(0)\ 



where the ground state of Ho is I'I'n ) with the 



r(0) 



,(0)| 



,(0)\ 



,(0)\ 



energy Eil" = (^^Wol*" > and |*^0 represents 
an excited state with the energy Em , whereas p 
denotes the number of lines connecting two given 
consecutive dots. 



„(0) 



Jl) 



„(2) 



(2) 



„(0) 



v(l) 



■ + zU- 



■ + Z(2-1)- 



TABLE I: Diagrammatic arrow line expressions of the non- 
vanishing elements cc^^^, including their multiplicities for a 
quadratic superlattice up to second hopping order. The coor- 
dinate number z is 2d for a d dimensional hypercubic lattice. 



(vl/(°'| 



and |*i°)) 



• ^WjV'I and iw;;'^ arc denoted by left-external and 
right-external lines, respectively, so in the diagram- 
matic representation of En there are some graphs 
which consist of s disconnected parts. The weight 
of these graphs has to be multiplied by the sign 

With these rules, we obtain within the line-dot represen- 
tation the perturbative expansion 



In our concrete case of the quadratic superlattice with the 
Hamiltonian in Eq. ([1]), the unperturbed Haniiltonian is 
given by 



ffo 



—UA (riA - 1) + —riB (riB ~ 1) 



-{^i + Afi)nA 



2 

flflB 



(A3) 



yielding the unperturbed energies ([25)) . whereas both the 
hopping and current terms are treated as a perturbation. 
Thus this leads to the arrow diagrams within the coef- 
ficient ckj"// , which can now be represented in terms of 

respective line-dot diagrams. Note that each term a^^, 
consists of exactly one creation operator (associated with 
Ji), one annihilation operator (associated with J*), and n 
hopping operators (associated with i"). For each arrow- 
line diagram we have to draw all possible topologically 
different line-dot diagrams. The sum of all these line- 
dot diagrams then gives the corresponding result. For 
example, the equation 



(A4) 



En = Ei"'> + A 




where i inside the dot stands for a particular sublattice, 
expresses exemplary how to transfer an arrow-line dia- 
gram into its line-dot representation. Following these 

(n) 

steps one obtains for the respective coefficients 0:211, ^^^ 
results (l30l)-(|32]), where the abbreviations ([281), (I2H1) are 
used. 
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